Cooling a charged mechanical resonator with time-dependent bias gate voltages 
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We show the possibility to cool a charged mechanical resonator (MR) down to nearly the vibra- 
tional ground state by elaborately tuning bias gate voltages on the electrodes, which couple the MR 
by Coulomb interaction. The key step is the modification of time-dependent effective eigen-frequency 
of the MR based on the Lewis-Riesenfeld invariant. With respect to a relevant idea proposed pre- 
viously [Li et at, Phys. Rev. A 83, 043803 (2011)], our scheme is simpler, more practical and 
completely within the reach of current technology. 
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Micro- and nano-mechanical resonators (MRs) have 
attracted much research interest due to their combi- 
nation of both classical and quantum properties [l|, 
together with broad impact on fundamental researches 
and applications 0, such as ultra-sensitive measure- 
ments approaching quantum limit 0, H| , observation of 
continuous-variable entanglement with mesoscopic ob- 
jects [H, Q i quantum information processing [?J , and bi- 
ological sensing P. 

The prerequisite of the research in these aspects is to 
cool the MRs down to their ground states in order to 
suppress detrimental influence from thermal fluctuations. 
Up to now, there have been many proposals for cooling 
the MRs |9M20] in either optomechanical or electrome- 
chanical systems, in which the resolved sideband cooling 
method enables the vibrational ground-state cooling of 
MRs. So far, the ground-state cooling of MRs has been 
achieved experimentally (2ll - [24| . 

Very recently, a different approach to fast ground- 
state cooling of MRs has been proposed [l6[ with time- 
dependent optical driving in a three-mirror cavity op- 
tomechanical system, in which the effective frequency of 
the MR can be changed like an 'optical spring' and the 
MR can be cooled down to nearly its ground- state un- 
der the control of time-dependent external optical driving 
fields. But the scheme seems challenging in experiments 
due to the requirement for both adiabatic evolution un- 
der Born-Oppenheimer approximation and very strong 
optical powers. 

We present in this work an alternative method, which 
is simpler but more practical, for cooling a charged MR 
by two electrodes via Coulomb interaction. Under the 
control of bias gate voltages, the MR behaves as a single- 
mode harmonic oscillator with its effective frequency tun- 
able like an 'electrical sprin g'. Our method is somewhat 
similar to that given in Ref. [la ] but without involvement 
of Born-Oppcnhcimcr approximation. Moreover, in con- 
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FIG. 1: Schematic diagram of the system, where a charged 
MR with (positive) charge Qmr couples to two identical elec- 
trodes via Coulomb interaction. The two electrodes are dis- 
tant by 2d and x is the deviation of the MR from the equi- 
librium position. Co is the capacitance of gate, and U is the 
tunable bias gate voltage. 



trast to Ref. [161 ] . our idea is more feasible using current 
experimental techniques due to manipulation of bias gate 
voltages. 

As schematically shown in Fig.[T] we consider a system 
where a charged MR is placed in the middle of two identi- 
cal electrodes, and couples electrostatically to these two 
electrodes. Each electrode takes the charge Q — CqU 
with Co and U being the capacitance and voltage of the 
bias gate, respectively. The Coulomb force between the 
MR and each electrode varies with the bias gate voltage, 
and the Coulomb potential is written by 
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where d is the equilibrium separation between the MR 
center-of-mass position and the electrodes with x the 
deviation of the MR from the equilibrium position, k 
(= l/47T£o) is the Coulomb constant in vacuum with Eq 
being the vacuum dielectric constant, and Qmr is the 
charge of the MR. 

After defining U = Uof(t) with the dimensionless time- 
dependent factor f(t) satisfying |/(f)| < 1, we may re- 
duce Eq. ([1]), under the condition of d ^> x, to 
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where the term 2kCoUoQMRf(t)/d has been ignored 
since it commutes with both the position and the mo- 
mentum operators of the MR, and makes no contribution 
to following discussions. 

Supposing the charged MR in the absence of the elec- 
trodes takes the Hamiltonian 
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(3) 



where x and p are the position and momentum operators 
of the MR with the bare eigen-frequency w m and the ef- 
fective mass m, we have the modified Hamiltonian under 
the control of the bias gate as 



H = H m + V c 

= i p2 + \ m ^ x2 



(4) 



where the effective frequency is w c ff = w m yl + r]f(t) 
with 77 = AkCoUoQMii/lmuj^d 3 ). 

Generally speaking, the dynamics of a harmonic oscil- 
lator (e.g. MR) governed by Eq. (@| can be obtained by 
the method of Lewis- Riesenfeld invariants p5l - l27| . Par- 
ticularly, for some special intermediate trajectories, the 
instantaneous quantum state of the harmonic oscillator 
at initial time instant can be the same as that at the fi- 
nal time instant, although the effective eigen- frequency 
of the harmonic oscillator has been changed significantly 
[27} . That is to say, we may have the final mean phonon 
number of the harmonic oscillator to be the same as the 
initial one, but with significant change of the effective 
temperature. This idea holds for arbitrary initial state, 
and has been used to cool the spatial motion of atoms 

25 1 and to cool the MR in cavity optomechanical system 

16 j. 

Our present idea is to achieve cooling of the MR by 
modifying time-dependent eigen-frequency via tuning the 
bias gate voltages on external electrodes. We note that 
the effective frequency uj e g of the MR can be easily 
changed by tuning the bias gate voltages U. In fact, w e fiE 
can be much smaller than u> m under negative U . More- 
over, u; ff can be even imaginary in order to accelerate 
the cooling process [l|| [25| . 

Specifically, considering some experimentally achiev- 
able parameters in charged MR systems j28l431} . such as 
x 10 9 N-m 2 /C 2 , tu m = 2ir x 134 kHz, m = 40 
3.15 ,um, U = 7.00 V, Qmr = \e\o- M R x s, 
C = 27.5 nF, and s = 0.08 
10 7 , meeting the approximate 



k = 8. 
Pg, d 

u mr = 1-25 x 10 13 /cm 2 , 
/1m 2 , we have r\ ~ 1.25 x 
condition in Eq. $2$. 

Similar to Ref. [l6j], our protocol for cooling MR 
comprises following two steps. 



Step I: Decreasing the mean phonon number with 
the increase of the MR effective frequency 

In the absence of bias gate voltages (i.e., f(t a ) = 0), 
we assume the MR with eigen-frequency u> e ff(t ) = ui m 



initially in a thermal equilibrium state 

H(tg) H{t a ) 

p(t a ) = e "b t /Tr(e ~ T ^ r ) 

at temperature T = 20 mK. So the corresponding mean 
thermal phonon number is given by 

n(t ) = 1/ {cxp[huj cS (t )/k B T} - 1} ~ 3100. 

With an arbitrary trajectory of oj e g(t) [or f(t)}, the 
MR is assumed to be in a new thermal state 
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at a later time ij with f(U) = 1, where the effective 
frequency is w e ff(Ai) — 3500oj m and the mean thermal 
phonon number is 



n(ti) 



1 



exp[/iw eff (t i )/(fc B T)] - 1 



0.47. 



(5) 



It is clear that the temperature of the MR at the 
time ti remains unchanged compared with the initial 
one, while the mean thermal phonon number of the MR 
is much less than the one at the initial time t Q due to 
enlargement of the effective frequency w e ff- The fact 
n(ti) ~ 0.47 < 1 means that the MR with enlarged 
effective frequency has been cooled down to nearly the 
vibrational ground state. 

Step II: Decreasing the MR effective frequency by 
keeping the low mean phonon number 

This step is accomplished via a special trajectory 
of f(t) evolving to the final time tf under the control 
of bias gate voltages satisfying two conditions: i) The 
final effective frequency is equal to the bare one (i.e., 
<jj e ff(tf) = uj m ), implying that the bias gate voltages are 
absent again at the final time tf (i.e., f(tf) = 0); ii) 
During the cooling process, the mean phonon number 
remains unchanged compared with that at time ti (i.e., 
n(tf) — fi(ti) ~ 0.47). Consequently, we have the MR 
with the bare frequency cooled to nearly the vibrational 
ground state. 

The solution to such a special trajectory of f{t) be- 
tween ti and tf governed by Eq. (Q| follows the Lewis- 
Riesenfeld invariant j25Tj27| . Using the inverse-invariant 
method [13, [HJ in Appendix, we obtain the control 
function of the trajectory of f{t) by tracing the Lewis- 
Ricscnfcld invariant, 
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where the dimensionless function b(t) is defined in Ap- 
pendix. 

The trajectory of f(t), starting from f(ti = 0) = 1 
and ending with f{tf) = 0, corresponds to the effective 
frequency of the MR decreasing from the end frequency 
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FIG. 2: (Color online) (a) Time evolution of f(t); (b) Time evolution of the corresponding effective frequency w e fi; (c) Time 
evolution of the effective temperature of the MR. In each panel, we consider three cases: tf — 0.5/w m (red solid line), 
tf — 1.0/wm (black dotted-dashed line), and tf = 2.0/oj m (gray dotted line). The parameter values are taken as ui m — 2n x 134 
kHz; m = 40 pg; a M R = 1.25 x 10 13 /cm 2 , Co = 27.52 nF, d = 3.15 (jtm, and s = 0.04/im 2 from [H-tH. 




FIG. 3: (Color online) (a) Time evolution of f(t) with 10% error(i.e., 0.9/(t) = (1 - 10%) f(t) and l.l/(t) = (1 + 10%)/(f)); 
(b)Time evolution of the effective frequency oi e // under the different trajectory of f(t); (c) Time evolution of the effective 
temperature of the MR under the different trajectory of f(t). In each plot, we consider two kinds: 0.0 f(t) = (1 — 10%)/(t) 
(thick line), and 0.9f(t) = (1 — 10%)/(t) (thin line). Except f(t), other parameters are the same as in Fig[2] 



in the first step (oj e s(ti) = lo ~ 3500u; m ) to the final 
frequency u e s(tf) — u m . We plot in Fig. [2] the con- 
trol function f(t), the corresponding instantaneous eigen- 
frequency ui e s(t) and the effective temperature T e ff. 

Using Eq. (5), we know that our cooling enables the 
MR from the initial temperature T c &(ti) = T c &(t ) = 
T = 20 mK down to the final temperature T e g(tf) = 6 
/xK, as plot in Fig.[2jc). Here the effective temperature is 
defined through n(t) = l/{exp[fuv e ff(t)/kBT e g(t)] — 1}. 

To estimate the imperfection in the trajectory of f(t), 
which corresponds to the experimental error in tunable 
bias gate voltages on electrodes, we give the trajectory of 
f(t) a ten percent fluctuation, i.e., (1 ± 10%)/(t), as in 
Fig. [UJa), and simulate in Fig. EJb) and (c) the changes 
of the instantaneous eigen- frequency uj c s(t) of the MR 
and the corresponding effective temperature of the MR. 
We find that, although the experimental deviation is as 
large as 10%, the effective eigen-frequency can still reach 
1.23w m for +10%/(t) (0.84w m for -10%/(t)) with the 
corresponding effective temperature 7 fiK for +10%/(t) 
(5 fiK for — 10%/(t)). In other word, in such cases, the 
MR can still be cooled down to nearly its vibrational 
ground state. 



We would like to point out that our scheme is only 
for an instantaneous cooling, which is different from the 
usual sideband cooling scheme achievable in steady state. 
This is because the MR is decoupled from electrodes at 
the end of our cooling process. So the MR can not be 
kept in low temperature for long time but be heated again 
by environment. 

In summary, we have proposed a practical protocol for 
cooling of a MR near to its vibrational ground state by 
controlling bias gate voltages on two nearby electrodes. 
Using the achievable experimental parameters and con- 
sidering experimental imperfection, we have shown the 
feasibility of the proposal by currently available technol- 
ogy. We argue that it is an alternative and practical idea 
for the MR cooling. 
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the National Natural Science Foundation of China (Grant 
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partly supported by the Research Funds of Rcnmin Uni- 
versity of China (10XNL016). 
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Appendix 

We present below how to obtain the special trajectory 
of f(t) between t = ti(= 0) and t = tf by means of 
the Lewis-Rcsisenfeld invariant. Based on the inverse- 
invariant method f2jl [32l . Eq. (J4j) satisfies the time- 
dependent invariant [2 



i T ^ 



2 l b(ty 



m 



(A.7) 



where loq = u) e g(ti) for simplicity, A = b{t)p — mb(t)x 
takes the role of the momentum conjugate to x/b(t), and 
the dimensionlcss real function b(t) follows 



b(t)+u 2 eS (t)b(t) 



j4_ 
b{tf 



(A.8) 



According to the boundary conditions at time t, and the 
finial time t f , we have 



b(U) = l, b(U) = 0, b(U) = 0, 
b(t f ) = x , b(t f ) = 0, b(t f ) = 0, 



(A.9) 



with x = a/wo / w e s (tf). So the simplest polynomial b(t) 
between tj (= 0) and tf is given by 

b(t) = 6(x - l)s 5 - 15(x - l)s 4 + 10(X - l)s 3 + 1 (A.10) 
with s = t/tf. 
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